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I Abstract 

i-C ' An analysis of arrival directions of extensive air showers generated by cosmic rays in the PeV energy 

I range and registered with the EAS MSU and EAS-1000 Prototype arrays reveals a considerable 

number of regions of excessive flux of cosmic rays. We present results of comparative analysis of 
regions found in the two data sets, estimate probabilities of their appearance, and discuss correlation 
of their locations with coordinates of possible astrophysical sources of PeV cosmic rays. 
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' 1. Method of The Investigation 

, The problem of the origin of cosmic rays (CRs) in the PeV energy range remains open for many years now. 
Attempts to find coincidences between arrival directions of CRs and their possible astrophysical sources is 
one of the routes of investigations on the subject. In the paper, we continue the earlier investigation [1] and 
present some of the results of a comparative analysis of arrival directions of extensive air showers (EAS) 
registered with the EAS MSU array in 1984-1990 and with the EAS~1000 Prototype ("PRO-1000") array 
^ in 1997-1999. 

! We selected 513 602 showers of the EAS MSU data set and 1 342 340 showers of the PRO-1000 data 

set for the purposes of the investigation. Showers registered with the two arrays have different number of 
charged particles Nc in a typical event. For the EAS MSU array, the median value of Nc is of the order of 
> ■ 1.6 X 10^ while that for the PRO-1000 array equals 3.7 x 10^. AU selected EAS have zenith angles < 45.7°. 
, The investigation is based on the method of Alexandreas et al. [2], which has been used by different 

research groups for the analysis of arrival directions of CRs. The idea of the method is the following. 
Every shower in the experimental data set obtains arrival time of another shower randomly. After this, new 
equatorial coordinates {a, 6) are calculated for the "mixed" data set thus providing a "mixed" map of arrival 
directions. The mixed map has the same distribution in S as the original map. Both maps are divided into 
' sufficiently small "basic" cells thus providing a way for comparison. 
, Mixed maps must be created multiple times in order to minimize the dependence of the result on the 

pseudo-random choice of arrival times. Finally, one calculates the mean of the mixed maps thus obtaining 
a "background" map. The underlying idea of the method is that this background map has most of the 
' properties of an isotropic background, and presents the distribution of arrival directions of cosmic rays that 
would be registered with the array in case there is no anisotropy. Thus, deviations of the real map from the 
background one may be assigned to a kind of anisotropy of arrival directions of EAS registered at the array. 

In our earlier analysis of arrival directions of EAS registered with PRO-1000, basic cells had the size 
1° X 1°, and the number of cycles of "mixing" varied from 100 to 200 [1]. In the present investigation, we 
reduced the size of basic cells down to 0.5° x 0.5° in order to improve accuracy of locating the regions of 
excessive flux (REFs) . This led to a slower rate of convergence of mixing. We thus had to perform a greater 
number of cycles of mixing. The difference in the number of EAS in basic cells of two consecutive maps is 
< 0.02 after 800 cycles for the EAS MSU data set, and after approximately 900 cycles for the PRO-1000 
data set. Thus we have chosen to run 1000 cycles of mixing for both data sets. 

Regions of excessive flux (REFs) of CRs were searched for in the following way. Adjacent basic strips 
in S (each 0.5° wide) were joined into strips with width AS = 3° . . . 30° with step equal to 0.5°. Each wide 
strip was then divided into adjacent cells of width Aa. After this, we calculated the number of EAS inside 
each of these cells for both experimental (iVrcai) and "background" (iVbg) maps. For each pair of cells, we 
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then calculated the value of significance S = {Nreai — A^bg)/\/-^bg- A cell was considered to be a region of 
excessive flux if 5 > 3. 

Rather unexpectedly, maps obtained after 1000 cycles of mixing with different random seeds in different 
runs lead to selection of slightly different sets of REFs. Due to this, we performed three 1000-cycle runs for 
the EAS MSU data set and two 1000-cycle runs for the PRO-1000 data set. Only those REFs found in all 
runs were selected for the following analysis. 

The method of Alexandreas et al. [2] does not put any requirements on the shape or size of regions in 
experimental and background maps to be compared. In what follows, we only consider so called "regular" 
regions. These are regions that have approximately the same area at different S and fixed AS, and have the 
shape of a square at 5 = 0°. To accomplish this, we used the standard formula Aa = AS/ cos S, where S is 
the mean value of S for the current strip, and Aa is rounded to the nearest half-integer number. 

The method of Alexandreas et al. does not provide a direct answer to the question about the probability 
of appearance of a REF. To solve the problem, we introduced a simple model based on the binomial distribu- 
tion. In the model, the number of trials equals the number of showers N in the data set under consideration, 
and an estimation of success (for a fixed region) equals p = N^g/N, where A'^bg is the number of showers in 
the region of the background map. The assumption is based on the fact that Ni,g may be considered as an 
expected number of showers in the region. Obviously, the probability of finding exactly Areai EAS in the 
region equals 

where v is a random variable equal to the number of successes in the binomial model, and C{N, N^cai) is the 
corresponding binomial coefficient. A more interesting probability to be considered is Pb = P{v < Aj-eai)- 

The condition 5 > 3 does not imply that it holds for sub-regions of a given region. Due to this, we have 
employed a number of additional quantities. One of the most usehil of them is a probability calculated from 
the following binomial model, which characterizes the degree of homogeneity of excess of the experimental 
flux in a given region over the background one. Let ^ be a random variable equal to the number of basic 
cells of the given region with an excess of EAS (over the background flux) . The number of trials equals the 
number of basic cells in the region (> 6 x 6 = 36). It is natural to assume the probability of "success" to be 
equal to 1/2. A probability to consider is P"*" = < N'^), where A''+ is the number of basic cells in the 
region with an excess of EAS in the experimental map over the background one. 

2. The Main Results 

Approximately 900 regular regions such that 5 > 3, P+ > 0.95 were found in the PRO 1000 data set, see 
Fig. 1. For all regions, max A5 = 17.5°, and the overwhelming majority of them have max A J < 12°. Nearly 
400 regions that satisfy the above two conditions were found in the EAS MSU data set, see Fig. 1. For all 
of them, AS < 10°. The probability Pb > 0.99812 for all REFs in both data sets. The number of EAS 
inside REFs differs from a few dozens up to approximately four thousand in the EAS MSU data set and 48.5 
thousand in the PRO-1000 data set. 

It is clearly seen that REFs found in the EAS MSU data set are, as a rule, more compact than that in 
the PRO-1000 data set. This may be related to the better accuracy in determination of arrival directions 
with the EAS MSU array. It is also clear that the number of coincidences or "embeddings" of REFs found 
in two data sets is comparatively small. In our opinion, the most interesting of them are the following: 

• a number of intersecting regions near the Supergalactic plane around a = 120°, S = 70° and a = 
183°, (5 = 30°; 

• a number of intersections of REFs located near the Galactic plane at a 323°-335°, and a ~ 350°-360°; 

• an embedding of a REF from the EAS MSU data set in a huge REF of the PRO-1000 set in the vicinity 

of a = 17°, S = 53°; 

• the only embedding of a REF from the PRO-1000 data set in a much larger EAS MSU region in the 

vicinity of a = 37°, 5 = 38°. 

Intersections of REFs located around a = 215°, (5 = 73° and a = 285°, 5 = 58° far from both Galactic 
and Supergalactic plans seem to be remarkable, too. One can also point out a number of pairs of REFs 

located close to each other. On the other hand, many REFs do not have "neighbours" in the other data set. 
In particular, there are no regular REFs satisfying the above criteria in the PRO-1000 set in the vicinity of 
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Figure 1: Regular regions of excessive flux such that S > 3 and P+ > 0.95 found in the PRO-1000 
(blue) and EAS MSU (red) data sets. For the sake of simplicity, joint boundaries instead of all 
individual REFs are shown in a number of cases. The arcs at the left and the right sides of the 
figure show the Galactic plane. The fl-like curve shows the Supergalactic plane. 

the Galactic plane in the direction close to the Galactic center (the lower part of the arc at the right-hand 
side of Fig. 1). 

Let us consider a relation between location of the REFs and coordinates of two traditional groups of 
possible Galactic sources of cosmic rays in the PeV energy range, namely supernova remnants (SNRs) and 
pulsars, see Fig. 2. 

Coordinates of the following nine SNRs are located inside or in the near vicinity of REFs from the EAS 
MSU data set (from top to bottom along the right-hand arc of the Galactic plane): G106.3-I-2.7, Cassiopeia A 
(3C461), CTB 104A, 3C434.1, 7-Cygni, G76.9-H1.0, CTB 87, G73.9-f0.9 and G69.7-hl.O. Almost all of them 
are of the shell type. For a considerable part of REFs in the EAS MSU set, there are pulsars with close 
coordinates and located at distances d = 0.33 ... 3 kpc from the Solar system. Among them, one can 
specify J1840-I-5640 {d = 1.70 kpc), J1838-hl650 (d = 1.84 kpc), J1709-I-2313 {d = 1.83 kpc), J1012-h5307 
{d = 0.52 kpc), J1115-h5030 {d = 0.54 kpc), J0814-t-7429 (d = 0.33-0.43 kpc), J0332-h5434 {d = 1.06- 
1.44 kpc) and a few other. At the same time, a few groups of REFs in the EAS MSU data set (including 
some of those with intersections with REFs in the PRO-1000 set) have neither SNRs, nor pulsars in their 
vicinity. 

As we have already reported earlier [1], the situation for the REFs in the PRO-1000 data set is qualita- 
tively the same. Namely, there are a whole number of Galactic SNRs (e.g., Cassiopeia A, the Crab Nebula 
(SN 1054), S147, 3C157, Tycho's Supernova (SN 1572), SN 1181 and G132.7-hl.3) and pulsars in the vicinity 
of REFs but still a number of regions do not have any of the above two classes of sources nearby. One of such 
REFs is a region located within a = 64.5° . . . 72.0°, 6 = 16° . . . 24.5°. It is especially remarkable because its 
boundaries are pretty close to the boundaries of an unidentified REF of 10-TeV cosmic rays found with the 
Milagro array [5]. To the best of our knowledge, the origin of this REF has not been explained yet [6]. 

3. Conclusions 

An analysis of arrival directions of EASs generated by cosmic rays in the PeV energy range and registered 
with the EAS MSU and EAS-1000 Prototype arrays reveals a considerable number of regions of excessive 
flux. Many of them have Galactic supernova remnants and pulsars located nearby but a number of REFs 
have none of the two types of cosmic ray sources in their vicinity. 

Only free, open source software was used for the investigation. In particular, all calculations were 
performed with GNU Octave [7] running in Slackware Linux. The investigation was partially supported by 
Federal contract No. 02.518.11.7073 and the Russian Foundation for Fundamental Research grant No. 08- 
02-00540. 
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Figure 2: Regular REFs found in the EAS MSU (upper figure) and PRO-1000 data sets and 
coordinates of Galactic SNRs (•) [3] and pulsars (o for pulsars with d < 3 kpc, A for those with 
3 < d < 6 kpc, X for the other) [4]. 
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